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Abstract 

The charge stabihty diagram for two coupled quantum dots containing up to two electrons is 
computed in magnetic fields. One- and two-particle Schrodinger equations are solved by exact 
diagonalization to obtain the chemical potentials and exchange energy in these systems. By an- 
alyzing the chemical potentials variation with external biases and magnetic fields, it is possible 
to distinguish between the weak and strong inter-dot couplings. The variation of the chemical 
potential curvatures and the double-triple point separations in the stability diagrams confirms the 
inter-dot coupling decrease with increasing magnetic fields. The computed exchange energies are 
also found to be significantly smaller than the values estimated from the stability diagram. 

PACS numbers: 73.21.La, 73.21.-b 
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Coupled quantum dots (QDs) are of particular importance for spin-based quantum com- 
putation because universal quantum logical gates (such as a Control-NOT gate) can be 
realized via the interaction between two quantum bits (qubits), i.e., the spins of two elec- 
trons, each trapped in one quantum dot.^ In such devices, the interaction between the two 
spins is proportional to the exchange energy J, which is equivalent to the splitting between 
the lowest singlet and triplet two-electron states. 

While extensive theoretical work focuses on the dependence of J on the system para- 
menters such as the inter-dot separation, the tunneling barrier between the QDs, and the 
external magnetic field,^^*^^ the charge stability diagram of coupled QDs^ has been studied 
to a lesser extent. Meanwhile, recent advances in experimental techniques have made it 
possible to study coupled QDs in the few-electron regime when each QD contains only one 
conduction electron (see, e.g., Refs. In this case the stability diagram becomes a 

powerful tool to study inter-dot coupling and electronic transport through double QD sys- 
tems. Analysis of the stability diagram and its evolution in magnetic fields allows one to 
estimate the values of the exchange energy as was demonstrated recently in the case of the 
two laterally coupled vertical QDs.— 

In general, in the stability diagram the boundaries between distinct stable charge states, 
i.e., between the states with fixed number of electrons A^^i and N2 in each of the coupled 
dots, are represented as functions of the two controlling gate biases, one for each dot.-^ These 
equilibrium charges are determined from the condition that the chemical potential of the 
QD structure iJi{Ni + N2) defined ass^ 



where Eg{N) is the ground state energy of the A^-electron state, is less than that of the 
leads (source and drain). 

In this paper, we numerically compute the stability diagram in coupled QDs with A^^i + 
N2 < 2 electrons in external magnetic fields, and investigate its properties for different 
inter-dot coupling strengths. The Hamiltonian for the coupled system is given by 



/i(iVi + N2) = EciNi + N2) - EciNi + N2- 1), 



(1) 



^^(ri,r2) = Horb + Hz, 



(2) 



Horb = h{ri) + h{r2) + C(ri, ra) 



(3) 
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2m* c 



(4) 



C(ri,r2) = eVe|ri-r2| (5) 
Hz = g^iBJ2B■Si (6) 

i 

Here, m* = 0.067me is the electron effective mass, e = 13.1 is the dielectric constant, 
g = —0.44 is the g-factor in GaAs, fiB is the Bohr magneton and A = ^[—By,Bx,0] is 
the vector potential for the constant magnetic field B oriented perpendicular to the QD 
plane [xy-plsme). The Zeeman effect is included in all our calculations except otherwise 
mentioned. 

The confinement in double QDs is simulated by the following model potential-^ 

where Vl and Vr are the depth of the left and right dots (equivalent to the QD gate volt- 
ages in experimental structuresS""^) which can be independently varied, d is the inter-dot 
separation, and R is the radius of the dot {R = 30 nm). 

The single-particle Hamiltonian Eq. Q and the orbital part of the two-electron Hamil- 
tonian Eq. (jH)) are diagonalized numerically to obtain the ground state eigenvalues which 
are used to compute the chemical potentials /i(l) and /i'^(2) [/i"^(2)] according to Eq. ((H) 
[note that Eg{0) = 0] and the exchange energy J as 

J = El{2) - Eg(2). (8) 

Above, the symbols S and T pertain to the singlet state and the lowest triplet state, respec- 
tively. 

We consider two coupled QD systems with different inter-dot separations li = 50 nm 
(strong coupling) and = 60 nm (weak coupling). The confinement potential for the two 
cases is plotted in Fig. ^a) where it is seen that aside from the increased lateral separation 
the barrier is also higher in the latter case (1.98 vs. 7.10 meV). Fig. ^b) and (c) show 
the chemical potentials /i(l) and yU'^(2) which in the absence of the magnetic field decrease 
almost linearly as Vl and Vn simultaneously increase. The quasi-linear decrease of the 
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chemical potentials with slightly larger slopes than shown was found for various magnetic 
fields as well.-^ 

We first analyze the dependences of the total energies Eg{1) and Eq{2) for one and two 
electrons on Vl and Vr shown in Fig. |21 We see that in the one-electron case (left column), 
the curvature of Eg{1) in the region where Vl and Vr are near each other {Vl ~ Vr) is 
larger for d = 60 nm than for d = 50 nm because of the weaker coupling between the QDs in 
the former case. We also note that for both values of the inter-dot distance, the curvatures 
of the energy contour plots increase along the main diagonal since the coupling between the 
dots decreases as Vl = Vr gets larger. 

However, when the two electrons populate the QD system, the situation becomes radically 
different: in the weak coupling case {d = 60 nm, bottom right), the total energy of the two 
electron system in the Vl ~ Vr region is almost linearly dependent on Vl (Vr), i.e., the 
curvature is vanishingly small, while for d = 50 nm (top right) the energy curves clearly 
exhibit a non-linear behavior with non-zero curvatures. The large overlap between the 
electrons in the strong inter-dot coupling case {d = 50 nm) is responsible for the smooth 
non- linear dependence of the energy on Vl (Vr). However, in the weak-coupling case {d = 60 
nm), the two electrons are well localized in the individual QDs by Coulomb repulsion and 
the large barrier between the dots, so that the potential change in one QD caused by the 
variation of Vl (Vr) does not affect the electron charge distribution but only acts as a 
constant addition to the total energy. This leads to a linear dependence of the total energy 
on Vl (Vr).— When the difference between Vl and Vr becomes sufficiently large to overcome 
the Coulomb repulsion, the two electrons move into one QD. This is accompanied by a change 
in the slope of the energy curves which become either horizontal or vertical as Vl (or Vr) 
no longer affects the total energy and which corresponds to the (0,2)/(2,0) regions on the 
stability diagram (not shown) .^^ We emphasize that the observed quasi- linear behavior of the 
total energy Eg{2) when Vl ^ Vr in the weak coupling regime {d = 60 nm) is physically 
different from the situation in two coalesced dots where both Eg{1) and Eg{2) are also 
straight lines perpendicular to the main diagonal in the Vl — Vr plane.® This is because in 
that case one deals with a single QD and changing Vl (Vr) modifies the total energy of the 
system. 

Fig. El displays the contour plots of (lower branches) and /i'^(2) (upper branches) 
as functions of Vl and Vr at zero magnetic field. We choose constant values of /^(l) = 
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fi^{2) = —18 meV for d = 50 nm and fi{l) = fi^{2) = —16.5 meV for d = 60 nm as the 
reference values of the chemical potential in the source/drain of the QD device.'^ In Fig. 
El we can recognize four regions corresponding to four stable charge states with different 
numbers of electrons in each dot [the numbers in the parentheses in each region give the 
number of electrons in the (left, right) QD] separated by the chemical potential contours 
and the main diagonal Vl = Vr. At the turning point on each branch along the main 
diagonal, three stable charge states coincide in terms of the total energy of the system. The 
distance between the turning points is the so-called double-triple point (DTP) separation 
(also called the anti-crossing separation) From Fig. El we also observe that the DTP 
separation AVl = AV/j = 5.00 meV in the ci = 50 nm case is significantly larger than 
the corresponding value AVl = AVr = 2.93 meV in the d = 60 nm case. Furthermore, 
the curvature of the branches around the DTP is smaller for d = 50 nm than for d = 60 
nm. According to the "classical" theory,® a smaller DTP separation (or equivalently a larger 
curvature of the chemical potential contour lines) indicates a weaker inter-dot coupling which 
is consistent with our findings. 

From the data in Table U we note that for d = 50 nm, the curvature (magnitude) /t(2) of 
the fi{2) curve is smaller than the curvature (magnitude) k(1) for /i(l), while in the d = 60 
nm case /t(l) < k(2). This peculiar behavior can be clarified by noting that both k(1) and 
/t(2) are determined by the differences between the corresponding curvatures of the total 
energy whose behavior in the voltage plane is discussed above. This indicates that in general, 
all being equal, in the weak-coupling regime the curvature of the chemical potential for two 
electrons is larger than that one for one electron, k(2) > while in the strong coupling 
regime, the opposite relationship k(2) < holds. 

In the presence of the magnetic field, the curvatures of the chemical potential contours also 
increase as can be seen in Fig. IHa) and (b) where we again plot the chemical potential con- 
tours for /i(l), At'^(2) and at constant reference values of /i(l) = /i'^(2) = fi^{2) = —18 
(—16.5) meV for d = 50 (60) nm at i? = 0, 3 and 6 T. Note the order of the contours for 
fi^{2) and /i^(2) at different magnetic fields. As the magnetic field increases, the contours 
shift from the lower left corner to the upper right corner because the single-particle eigenen- 
ergies increase.^ In addition to the curvature increase, the DTP separation becomes smaller 
at larger magnetic field for both singlet and the lowest triplet states [for the detailed expla- 
nation of this effect, see the discussion on Fig. El^a)]. From the changes in the curvature and 
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DTP separation, one concludes that the magnetic field indeed induces a quantum mechani- 
cal decoupling between the two dots and results in magnetic localization of electrons in each 
dot. By comparing the chemical potential curves in Fig. ^a) with the corresponding ones 
in Fig. Ilfb), we see that in the latter case the chemical potential contours have much larger 
curvatures than in the former case (see Table |l] for details) due to the increased inter-dot 
decoupling and for each value of the magnetic field the DTP separation for d = 60 nm is 
more than 60% smaller than for d = 50 nm. 

From Table m it is also shown that the curvatures k,{1) and k(2) progressively increase as 
the magnetic field becomes larger. This is due to enhanced localization of electrons caused 
by the magnetic field. The magnetic localization in the weak coupling case became prevalent 
at lower fields than in the strong coupling situation [see lower insets of Fig. Ej, which is 
manifested by a more rapid increase in the curvature of chemical potential contours. 

FiguresEfa) and (b) show the extracted DTP separation along Vl (or Vr, Vl = Vr) axis as 
a function of magnetic fields for d = 50 nm and 60 nm inter-dot separations, respectively. In 
each plot the data are shown for the singlet and lowest triplet states. Note that at i? = the 
DTP separation for the singlet state is smaller than that for the lowest triplet state because 
the singlet is the ground state, while at larger B fields, the lowest triplet state becomes the 
ground state and the order of the DTP separations is reversed. In both (a) and (b), the DTP 
separation for the lowest triplet state decreases faster with B fields than that for the singlet 
state. This is because the DTP separation is proportional to /i(2) — /i(l) = Eg{2) — 2Eg{1) 
for a fixed Vl = Vr on the main diagonal of the stability diagram (see Figs. ^ and El). For 
the singlet state, Eg{2) does not change with the B field while Eg{1) decreases with the B 
field due to the Zeeman effect, therefore the Zeeman contribution to /i(2) — fi{l) increases 
with the B field. For the triplet state, the Zeeman contributions to Eg{2) and 2Eg{1) cancel 
out, and /i(2) — /i(l) is not affected by the B field. The decrease of the DTP separation 
in the magnetic field was also recently observed experimentally.— The upper (lower) inset 
in each figure shows the corresponding exchange energy J as a function of the magnetic 
field calculated by Eq. (jH} with (without) the Zeeman effect. In both cases, the Zeeman 
effect induces a linear depenedence of J on B. However, in (a) given the strong coupling 
between the dots, the orbital contribution to J dominates at low B fields before being 
overcome by the Zeeman induced decrease at higher field; in (b), J is totally dominated by 
the Zeeman contribution, which decreases linearly with the B field. Comparison of the B 
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field dependences of tlie DTP separation and excfiange energy in tfie absence of the Zeeman 
effect in Fig. Elshows tliat tlie latter saturates at much lower values of the magnetic field than 
the former. This is because the DTP separation is determined by the Coulomb interaction 
between electrons which decreases as the electrons become localized by the magnetic field in 
individual dots (within the Heitler-London approximation, this decrease is proportional to 
B~^, Ref. Q), while the exchange energy in absence of the Zeeman effect approaches zero 
much faster than the Coulomb interaction since it is proportional to the overlap between 
the individual electron wave functions that decays exponentially fast in strong magnetic 

It is also interesting to compare the exact values of the exchange energy (see the insets 
in Fig. EI) with those extracted from the stability diagrams in magnetic fields using the 
Hubbard model.-^ According to this model, Jest = 4t^/ (Vintra — Vinter) where 2t is the tunnel 
(symmetric-asymetric) splitting, Vintra and Vinter are the intra-dot and inter-dot Coulomb 
interactions. From the data shown in Fig. we estimate the value of the inter-dot Coulomb 
interaction ^ 3.4 (2.0) meV for d = 50 (60) nm, which is given by the DTP separation 

(for the lowest triplet state) in the limit of large magnetic fields. These numbers are in good 
agreement with the corresponding expectation values (C(ri, r2)) of the Coulomb interaction 
matrix (3.5 and 2.2 meV, respectively) obtained from direct calculations, thereby confirming 
electron localization and QDs decoupling. Since at zero magnetic field, the DTP separation 
is equal to 2t + Vinter, we obtain 2t^^^^'^^ ^ 1.6 (0.7) meV which is consistent with the energy 
differences between the two lowest single-particle levels of 1.9 (0.4) meV. As Vintra ~ 8 
meV is given by the electron addition energy in one QD which is the distance between the 
" corners" of the linear region where single electron re-localization occurs from one dot to the 
other in the N = 2 energy diagram (see Fig. El), the estimated values of the exchange energy 
become Jest^^^ ~ 0.6 (0.08) meV. These numbers are of the same order as the numerically 
exact values of 0.24 (0.012) meV, but they both significantly overestimate the computed 
data, and therefore, can only be used as a general guideline to gauge the magnitude of 
the exchange coupling in double QDs. The overestimation is due to the difference between 
Coulomb energies in the singlet and triplet states that lowers the exchange energy,— but 
which is not taken into account in the simple Hubbard model. 

In summary, we computed the stability diagram for model double QD systems populated 
with up to two electrons in magnetic fields using numerically exact diagonalization of the one- 
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and two-electron Hamiltonian. Two inter-dot separations d — 50 and 60 nm corresponding 
to strong and weak inter-dot coupling were considered. We found that in the weak-coupling 
regime the curvature of the chemical potential /i.(2) is larger than that one of while 
in the strong-coupling case the situation is reversed. Hence, by analyzing the chemical 
potential variations caused by external biases and magnetic fields, it is possible to distinguish 
between strong and weak inter-dot coupling, even if the curvatures are of the same order 
in both cases. The evolution of the stabihty diagrams in magnetic fields conforms to the 
general idea of enhanced electron localization with increasing field strength. The exchange 
energies extracted from the stability diagrams showed that these values are significantly 
overestimated when compared with numerically exact data. 

This work is supported by DARPA QUIST program through ARO Grant DAAD 19-01- 
1-0659. The authors thank the Material Computational Center at the University of lUinois 
through NSF Grant DMR 99-76550. LXZ thanks the Computer Science and Engineering 
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FIG. 1: (Color online) (a) The confinement potential for d = 50 nm (red) and d = 60 nm (blue) 
at Vl = Vr = 25 meV. Chemical potentials ^(1) (solid) and fJ-^{2) (dashed) vs. Vl = Vr for (b) 
d = 50 nm and (c) d = 60 nm. In both of them the horizontal line indicates the values of the 
chemical potential at which the contours in Fig. |21are drawn. 
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FIG. 2: (Color online) Surface (contour) plots of the total energies for = 1 (left column) and 
N = 2, singlet (right column) at (i = 50 nm (top row) and d = 60 nm (bottom row). In all plots 
the energies decrease from the lower left to the upper right corner. 
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FIG. 3: (Color online) Contour plots of the chemical potentials and /U'^(2) as functions of 
Vl and Vr for d = 50 nm (red) and d = 60 nm (blue). The turning points on the contour lines 
are indicated by solid dots and the dotted line is a guide for the eyes along the main diagonal 
(^L = Vr)- The numbers on the left (right) within parentheses give the electron number in the 
left (right) dot [the (0,0) region is located at the lower left corner below the //(I) branch for d = 50 
nm]. 
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FIG. 4: (Color online) (a) Contour plots of the chemical potentials //(I) (single electron state, 
lower branches, solid curves), /^'^(2) (two-electron singlet state, upper branches, solid curves), and 
^^(2) (the lowest two-electron triplet state, upper branches, dashed curves) at different magnetic 
fields for d = 50 nm. (b) Same as (a) but at d = 60 nm. In the case of d = 60 nm and B = 0, the 
contour lines for ^'^(2) and /i"^(2) are indistinguishable on the scale of the figure. 
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FIG. 5: (Color online) Double-triplet point separation along Vl (or Vr, Vl = Vr) axis as a function 

of the magnetic field B for (a) d = 50 nm and (b) c? = 60 nm. The data for singlet and the lowest 

triplet states are labeled by "A" (blue curve) and "+" (red curve), respectively. The upper (lower) 

inset in each figure shows the exchange energy J as a function of the magnetic field with (without) 

the Zeeman effect. The data in the insets are obtained at Vl = Vr = 25 meV. 
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TABLES 



TABLE I: Comparison of the magnitude of the second order derivative k{1) and k{2) of the 
chemical potential curves and /x(2) near their turning points for various inter-dot separations 
and magnetic fields. k^{2) and k^{2) denote the values for the singlet and the lowest triplet state, 
respectively. 
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0.400 
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6.145 [7.463] 
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0.714 


0.631 [0.659] 


15.706 


27.230 [37.141] 



15 



